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Abstract. We study the Doppler-broadened absorption of a weak monochromatic 
probe beam in a thermal rubidium vapour cell on the D lines. A detailed model of the 
susceptibility is developed which takes into account the absolute linestrengths of the 
allowed electric dipole transitions and the motion of the atoms parallel to the probe 
beam. All transitions from both hyperfine levels of the ground term of both isotopes 
are incorporated. The absorption and refractive index as a function of frequency are 
expressed in terms of the complementary error function. The absolute absorption 
profiles are compared with experiment, and are found to be in excellent agreement 
provided a sufficiently weak probe beam with an intensity under one thousandth of the 
saturation intensity is used. The importance of hyperfine pumping for open transitions 
is discussed in the context of achieving the weak-probe limit. Theory and experiment 
show excellent agreement, with an rms error better than 0.2% for the D2 line at 16.5°C. 
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1. Introduction 

In this paper we develop a model that allows us to predict quantitatively the absorptive 
and dispersive properties of rubidium vapour probed in the vicinity of the D lines 0, 
and compare the predictions with an experimental study of the absolute Doppler- 
broadened absorption spectrum. In alkali-metal atoms the D lines have a high oscillator 
strength, and from an experimental perspective rubidium and cesium are ideal species 
as they have sufficient room temperature vapour pressure to yield large absorption 
(10-50%) in cells of modest length (5-7 cm). In addition, for these atoms the D 
lines occur at wavelengths where inexpensive and reliable diode laser sources exist. 
These transitions are frequently studied in atomic physics; the D2 transition is used 
extensively in laser-cooling experiments [HIS], whereas non-linear optical processes such 
as electromagnetically induced transparency [3] and chip-scale atomic magnetometers [1] 
are realized with the Di transition. 

Controlling the propagation of light through a medium by modifying its absorptive 
and dispersive properties is a flourishing area of research [5l[6l[7]. Having a model which 
calculates the absorption and refractive index of a Doppler-broadened medium is useful 
for predicting the magnitude of relevant pulse propagation properties, a topical example 
being the group delay which enables construction of an all-optical delay line [8], [9] or 
a slow-light Fourier transform interferometer [TOj. A model which predicts absolute 
absorption has a number of applications e.g. in analysing EIT spectra [HI [12]. Absolute 
absorption spectroscopy can yield the number density of the sample being studied and 
has many applications in physics, chemistry, metallurgy and industry [131 [E]; when 
applied to the measurement of different spectral lines the populations of occupied levels 
are revealed, from which a temperature can be deduced [I5]. In addition, Rb and 
Cs spectroscopy is frequently used for generating a signal used for frequency reference 
( "locking" ) of a laser. Understanding the evolution of the absorption profile aids in the 
decision of, e.g. which temperature to use to maximise the signal ^16j. The importance of 
using a weak probe in order to maximise the absorption will be highlighted below. Many 
laser-lock schemes have signals which have a non-trivial dependence of signal amplitude 
on absorption, e.g., polarization spectroscopy [17J or the dichroic atomic vapour laser 
lock (DAVLL) [18] . The choice of which probe power to use is a trade-off between two 
competing effects: a weak probe beam ensures that the largest absorption is obtained, 
whereas higher probe power gives a better signal-to-noise ratio. 

The aim of this work is to provide a detailed model of the absorption and refractive 
index for the Doppler-broadened Rb 55* — > 5P transition, and to compare the expected 
spectral dependence of the absorption with experimental data. The structure of the 
remainder of the paper is as follows: Section [2] explains how to calculate the expected 
Doppler-broadened spectra for Rb atoms on the D lines; Section [3] describes the 
experimental apparatus and details of the methodology used to measure the absorption 

I For an alkali metal atom the D2 transition is n'^Si/2 —>■ where n is the principal quantum 

number of the valence electron, and the Di transition is n^5i/2 n'^Pi/2- 
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profile; Section H] presents and discusses the results, and finally, in Section El we draw 
our conclusions. 

2. Calculating the absorption coefficient of a Doppler-broadened medium 

2.1. The atomic absorption cross-section 

The absorption of monochromatic light as it propagates along the direction through 
a uniform density atomic vapour is given by the Beer-Lambert law: 



where I{z) is the intensity of light at position z inside the medium with an absorption 
coefficient a(z/, T), and Jq is the beam intensity at the entrance of the medium. The 
absorption coefficient is dependent on the frequency, z/, of incident light and the 
temperature, T, of the medium. We assume that the probe beam is sufficiently weak 
that the absorption coefficient is independent of intensity. A full discussion of how weak 
the light has to be for this simplification to be valid is given in section 14.11 

The transmission, T, of a beam through a medium of length L is defined as 



^0 

In general, a medium consists of multiple species, each with multiple transitions. A 
beam of light will interact with all species according to 



where the total absorption coefficient, Sctj, is the sum over Oj, the absorption coefficient 
for each transition for each species. 

The macroscopic absorption coefficient of the medium can be written [19] in terms 
of a, the microscopic atomic absorption cross-section and Af, the number density of the 
atomic gas, a = Ma. There are two reasons why the medium's absorption coefficient is 
temperature dependent: (i) the atomic cross-section is influenced by the Doppler width, 
proportional to the square root of the temperature; and (ii) the number density is a 
strong function of temperature. Doppler broadening of the spectral lines is dealt with 
in section 12.51 and the temperature dependence of the atomic density in Appendix A. 

We label each hyperfine state of the atom with the usual angular momentum 
quantum numbers \Fg.,'mFg) for the ^6*1/2 term, and |Fe,mi?J for the ^^3/2 or ^Pi/2 
term, where the subscript g{e) denotes the ground (excited) state. For a multi- level 
atom, such as rubidium, the calculation of the atomic cross-section is in two parts: 
first, the relative linestrengths among the different \Fg.,mFg) \F(..,mp^) transitions 
are calculated, then the absolute value is deduced. These calculations are facilitated 
by initially assuming the atoms are at rest, with the manifestation of atomic motion 
(Doppler broadening) incorporated later. We neglect pressure broadening in these 
calculations. Gorris-Neveux et al. measured the Rb-Rb coUisional self-broadening 




(1) 




(2) 



T = exp[— (Saj)L] 



(3) 
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to be of the order of 10~^ Hz cm^ [20]; for the temperature range spanned in this 
work the pressure broadening is at least four orders of magnitude less than the natural 
broadening. For temperatures greater than approximately 120°C the self-broadening 
becomes comparable to the natural width. 



2.2. Transition frequencies 

In order to predict the absorption spectrum the relative spacing of the hyperfine- 
resolved energy levels for both Rb isotopes are needed. Zero detuning frequency 
for D2 (-Di) is set to be the centre of mass frequency of the 5s^S'i/2 — > 5p^P3/2 
(5s^S'i/2 — ^ 5p^Pi/2) transition in the absence of hyperfine splitting, taking into account 
the natural abundance of each isotope. The atomic energy level intervals were obtained 
for D2 from [211 [22l [23] and for Di from [24J . The positions of the atomic transitions 
relative to the centre of mass for D2 (384,230,426.6 MHz) and Di (377,107,407.299 MHz) 



are listed in tables 1 (a) and 1 (b) respectively. 



Table 1. Transition frequencies for the (a) D2 line, and (b) Di. 
(a) (b) 
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2. 3. Relative linestrength factors 

The strength of the interaction between an atom and near-resonant electromagnetic 
radiation is characterized by the dipole matrix elements. The dipole matrix element of 
the transition between states \Fg,mpg) and \Fe,mpe) is {Fg,mpg\erg\Fe,mpe) ■ In order 
to calculate this matrix element, it is possible to factor out the angular dependence and 
write the matrix element as a product of Wigner 3 — j and 6 — j symbols and a reduced 
matrix element [251 [26] . Thus, 

{Fg,mpg\er,\F,,mp,) = (_i) 2F.+x+j,+j.+L,+5+m^,+i^^^| |^^| |^^^ 

X ^i2Fg + l)(2Fe + 1)(2J, + 1)(2J, + 1)(2L, + 1) 



Absolute Rb absorption 



5 




[• (4) 

Here F,X, J, L, S, and mp are the angular momentum quantum numbers, and q is the 
integer change in during the transition. T, the nuclear spin, has the value | and | 
for ^^Rb and ^^Rb respectively. S, the electron spin, has the value |. The 3 — j symbol 
is the term contained in the large round brackets, and the Q — j in curly brackets. Note 
that the 3 — j symbol is non-zero for mpe = ^Fg + according to the usual definition 
of q. {Lg\\er\\Le) is the reduced matrix element, and can be expressed in terms of the 
wavelength of the transition. A, and the decay rate of the excited state, T. By calculating 
the Wigner coefficients and prefactors, equation (jl]) reduces to 

{Fg,mFg\erg\Fe,mFe) = c„^(Lg| |er| |Le) = c^pd, (5) 

where Cmp is a coefficient that determines the transition strength of a particular 
transition, and is dependent on the initial and final states of the transition. 

The strength of a transition is proportional to the square of the transition matrix 
element, thus the transition strength is c^^d?- Each hyperfine transition is degenerate 
in F (since we are assuming zero magnetic field). The total transition strength of 
the hyperfine transition Fg — > F^ is denoted by = Sc^^, where is the sum of 
transition strengths c^^ of each Zeeman transition in the hyperfine manifold. These Cf^ 
coefficients have been calculated for linearly polarised light (g = 0), and are tabulated 
Appendix B 



m 



2.4- Absolute absorption coefficient 

The reduced matrix element, d, can be calculated using the expression for the decay 
rate [19] 

(Jg||er||Je) can be written in terms of (Lg||er||Le) via the relation 

(J,||er||Je) = (-l)-^^+^«+^+^(L,||er||L,) (7) 

Lg 1 I 

Je Jg S ] 



X J(2Je + l)(2L, + l) 



The Wigner 6 — j coefficient and prefactor, both of which are independent of the F and 
mp quantum numbers, can be calculated for the D2 line. Thus 

{Jg = l/2||er||Je = 3/2) = ^^{Lg = 0||er||Le = 1). (8) 
Substituting ([8]) into ([6]) and rearranging. 



d={Lg = Q\\er\\U = l) = V3f-^^. (9) 
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For the Di line a similar analysis leads to the same result as equation The reduced 
dipole matrix element for the fine structure splitting should be identical for the D lines. 
However, we have used experimentally measured values for the wavelength and decay 
rates: A = 780.241 nm [27j and T = 27r x 6.065 MHz [28] for D2, and A = 794.979 nm [29] 
and r = 277 X 5.746 MHz [28] for Di. This yields d = 5.177 eao for D2, and d = 5.182 eoo 
for Di, where oq is the Bohr radius. 

2.5. Including atomic velocity 

The thermal velocity of atoms along the axis of the probe beam is given by the well- 
known Maxwell-Boltzmann distribution. It is Gaussian in nature, with a 1/e width of 
u = ^ 2kBT / M , where T is the temperature, and fee is the Boltzmann constant and M 
is the atomic mass. It is this longitudinal motion that leads to Doppler broadening of 
the absorption spectra. Let the angular frequency of the laser be ui^, and that of an 
atomic resonance be Uq. The angular detuning. A, is defined as A = — ujq. For an 
atom moving along the direction of propagation of the probe beam we incorporate the 
Doppler effect by simply replacing the detuning hj A — kv, where k is the magnitude of 
the wavevector of the light, and v is the atomic velocity. We assume that the experiment 
is conducted in the weak-probe limit, i.e. the laser intensity is sufficiently low that 
optical pumping processes which redistribute population amongst the hyperfine levels 
of the ground term do not occur during the transit of an atom across the finite beam 
width. The transverse motion of atoms can therefore be neglected. 



2.6. Electric Susceptibility 



The susceptibility, x{^): encapsulates both the absorptive and dispersive properties of 
a medium. For a medium composed of atoms at rest, the susceptibility for the transition 



Fe is given by 



XF,F.(A) = C^rfW— /r(A) 



(10) 



where Cpd"^ is the transition strength of the hyperfine transition, Af is the number 
density, and /r(A) is a lineshape factor derived from the steady state solution to the 
optical Bloch equations of a two-level atom, in the absence of Doppler-broadening. 

As the atomic dipoles are not in phase with the driving light field x(^) is, in 
general, a complex function; the real part characterises dispersion, and the imaginary 
the absorption. The susceptibility for atoms with velocity v along the beam propagation 
direction is given by [19] 



fr{A-kv) 



A-kv 



-1 



[A - kv) 



A — kv\2 
2 



+ r 
2 



(11) 

A — A;f \2i -1 
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/]? and /p denote the real and imaginary parts of /r respectively, with /j? having a 
characteristic dispersion profile, and // being the Lorentzian absorption profile expected 
for an homogeneously (natural) broadened system. 

The magnitude of the susceptibility depends on the strength of the transition in 
question, which is simply a prefactor. Hence we define for convenience s(A), which is 
directly proportional to but is independent of the specific atomic transition. By 

integrating over the atomic velocity distribution, one obtains the Doppler-broadened 
lineshape 



r+oo 

s{A) = I fr{A-kv)xg^{v)dv, (12) 



where 



9u[v) = exp 



(13) 



u 



is the normalised Gaussian, with 1/e width u. 

Making the substitutions y = A/ku, x = v/u and a = T/ku, equation f|T2|) becomes 

/ + 00 
fa{y -x) X g{x)dx, (14) 
-CXD 

where s{y) is in units of {ku)~^, and g{x) has been normalised, with a dimensionless 
width of 1. This is in the form of a convolution integral, and can be re- written as 
s{y) = fa{x) ® g{x). Separating the real and imaginary parts of /a, 

s{y) = fai^) ® 9{x) + ^flix) ® 9{x), (15) 

and hence 



s'^iy) = Re[/f (x) ® g{x) + ifi{x) ® g{x)] (16) 
= fa{x)®g{x) 

s\y) = Im[/«(x) ® g{x) + ifi{x) ® g{x)] (17) 
= fi{x)(S>g{x), 

where the fact that the convolution of two real functions is real has been used. The 
imaginary part is related to the absorption coefficient, and is the well-known Voigt 
function, being the convolution of a Lorentzian and a Gaussian function. 

Using the convolution theorem of Fourier transforms, the convolution of two 
functions can be re-written as 

S{y) = F,{x) X G{x), (18) 

where capitals denote the Fourier transform of a function, and the tilde denotes the 
reciprocal variable. The advantage of using this method is that the Fourier transform 
S{y) is simply the product of two functions which can be calculated analytically. s{y) 
can then be produced from S{y) by taking its inverse Fourier transform. 
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For the case of the Voigt profile, the Fourier transforms of (x) and g{x) are 



F^{x) = Trexp 



a , 

2 



\x\ 



G(x) = exp 

Taking the inverse Fourier transform of their product results in a Voigt profile of 



v^ei(a-2.)^^Erfc 



- W 



+ e^=^"^Erfc 



+ iy 



(19) 
(20) 

(21) 



where s^{y) is in units of (ku)"^, and Erfc[2;] denotes the complementary error function 
of z. 

The real part of the susceptibility is related to the refractive index, and can also 
be expressed in terms of the complementary error function: 



s«(^)^i^^(el(«-2,)2Erfc 



ty 

2 ^ 



a 

^ + ^y 



(22) 



s^{y) is in units of (ku)'^ and, despite the prefactor of i, is entirely real. s^{y) can be 
differentiated with respect to y to arbitrary powers to evaluate, e.g., the group refractive 
index. 



2.7. Absorption Coefficients 

The absorption coefficient can be obtained from the imaginary part of the susceptibility, 
X(A), via 

a{A) ^ klm[x{A)], (23) 

where k is the wave number of the probe beam. Im[x(A)] has the form of a Voigt 
profile, s^(A), multiplied by prefactors which depend on the properties of the resonant 
transition. 

The width of the Voigt profile is characterised by a single parameter, a: the ratio 
of the widths of the Lorentzian to the Gaussian profiles. The width of the Lorentzian, 
r, is the full-width at half-maximum (FWHM) of the hyperfine-free atomic transition, 
r is identical for all hyperfine transitions and Zeeman sublevels within the hyperfine- 
free manifold, and is also equal for different isotopes of an element. The width of the 
Gaussian profile is proportional to the width of the Maxwell-Boltzmann distribution, u, 
and is a function of temperature and isotopic mass. Gonsidering all of the above, the 
width of the absorption profile of every hyperfine transition for a particular isotope is 
identical. 

The height of the Voigt profile depends on two factors: the forms of the Lorentzian 
and Gaussian functions, which are identical for all transitions for a given isotope; and 
the transition strength of a particular transition. 

All transitions for a particular isotope can be represented by a single Voigt profile, 
which is then centred on the relevant transition frequencies, and multiplied by the 




Figure 1. Plots of the transmission througli a vapour cell of length 75 mm as a 
function of linear detuning, A/27r. Plot (a) shows D2 at 20°C, and (b) shows Di at 
30°C. The blue lines [i) show the transmission for the transitions ®^Rb Fg = 2 ^ F^, 
the magenta (ii) *^Rb Fg = 3 ^ Fg, the green {in) ^^Rb Fg = 2 ^ Fe, and the orange 
(iv) ^^Rb Fg = 1 ^ Fe- The solid lines show the transitions between hyperfine states 
Fg^ Fe^ Fg + 1, dashed Fg^ Fe = Fg, and dot-dash Fg^ Fe = Fg-l. The black 
line shows the total transmission through the cell. Zero detuning corresponds to the 
weighted centre of the line. 



relevant transition strengths. Hence, recalling equations ( ITOl) and (l23ll . the absorption 
profile for a particular hyperfine transition Fg — > F^. is 

.,,,.(A) = .I,n|,(A)l = kCU^M^^^^-?^. (24) 

Here, 2 (2X + 1) is the degeneracy of the ground state of the particular isotope (12 for 
^^Rb, 8 for ^"^Rb). The degeneracy appears as we are assuming that the population is 
evenly distributed amongst the ground state Zeeman sublevels (at room temperature 
the Boltzmann factor is 1 for the two different Fg hyperfine states, and reduces the 
population of excited states to a negligible level). The expected transmission profile for 
the vapour cell can then be calculated as a function of detuning. 

Figure [1] shows the predicted transmission spectrum for rubidium vapour in a 
75 mm-long cell for (a) D2 at 20°C, and (b) Di at 30°C. The contributions of the 
individual Fg F^ transitions are shown, in addition to their combined total. For 
both D lines the ground state hyperfine splitting is larger than the Doppler width of 
~ 0.5 GHz. For the D2 line, the excited hyperfine splitting of both isotopes is smaller 
than the Doppler width; consequently four composite lines are observed. For the Di 
line, the excited state splitting for ^^Rb is smaller than the Doppler width, whereas the 
splitting for ^^Rb is larger; hence six composite lines are seen. 




Figure 2. Schematic of the experimental apparatus. Light from an external cavity 
diode laser (ECDL) passes through an optical isolator (01) and impinges on a beam 
splitter (BS). A fraction of the beam passes through a Fabry-Perot etalon onto a 
photo detector (PD). Another fraction of the beam is expanded in a telescope and 
passes through aperture A. Mirrors (M) steer the beam and another beam splitter is 
used to make pump and probe beams which cross at a small angle in a Rb vapour cell. 
The probe beam is incident on a photo detector, and neutral density filters (ND) are 
used to give independent control over the pump and probe beam powers 

3. Experimental methods and results 
3.1. Experimental apparatus 

We now test the accuracy of the prediction experimentally. A schematic of the 
experiment is shown in figure [2l External cavity diode lasers were the source of light 
(Toptica DLIOO at 780.2 nm and 795.0 nm for D2 and Di respectively). A fraction of the 
output beam was used as a probe beam for rubidium vapour in a 7.5 cm cell. A portion 
of the light was also sent into a Fabry-Perot etalon. A telescope was used to expand the 
probe beam before the cell. Before the cell the beam had a radius of (2.00±0.05) mm. 
The cell could be heated to change the vapour pressure of rubidium and hence the 
opacity. A thermocouple was used to measure the approximate temperature of the 
cell. No attempt was made to null the laboratory magnetic field. A pump beam 
generated sub-Doppler spectral features to provide a frequency reference. The crossing 
angle between probe and counterpropagating pump within the vapour cell was 6 mrad. 
Neutral density filters were used to give independent control of the pump and probe 
powers. The Fabry-Perot etalon was used to assist with calibrating and linearizing the 
frequency scan. A plane-plane cavity was used, with a separation of the mirrors of 
25 cm, with a free-spectral range of 0.60 GHz. The probe beam was incident on a photo 
detector comprising a simple current-to-voltage circuit designed to output a voltage 
linearly proportional to the incident power. 
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3.2. Scaling the frequency axis. 

For the D2 line the frequency axis of the laser scans were linearized by use of the 
etalon transmission peaks. In order to generate atomic frequency markers on a scale 
narrower than the Doppler-broadened features pump-probe spectroscopy was employed. 
By counterpropagating a pump beam with the weak probe it is possible to generate 
sub-Doppler saturated-absorption and hyperfine-pumping spectra [30j. For each isotope 
ground state Fg, three transitions Fg ^ Fe = Fg, Fg±l are resolved, and so-called cross 
over resonances [30] are seen halfway between each resonance. In this way it is possible 
to obtain 24 atomic resonances. 
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Figure 3. (a) Deviation of the measured spectral line frequencies from their expected 
positions for the D2 line before linearization of the laser scan. The black crosses mark 
the measured positions of the sub-Doppler spectra, (b) Linear fitting of the expected 
to the measured frequency after linearization. The red line shows a linear relationship 
between the two axes, with a gradient of 1 and an intercept of zero. The inset shows 
the deviation of the measured spectral line frequencies from their expected positions. 



Figure [3](a) shows a plot of the difference between the measured and expected 
detunings of the 24 atomic resonances before linearization. The expected detunings were 
obtained from table 1(a) The relatively large deviations from zero are seen to have a 



polynomial relationship with expected frequency. Figure [3]^b) shows a plot of measured 
detuning of the atomic resonances versus the expected detuning after linearization. The 
solid line has slope 1 and passes through the origin. The inset shows the deviation 
between measured and expected frequency. It can be seen that each atomic resonance 
is within 5 MHz of this ideal fit over a span of 8 GHz. The residual deviations are a 
consequence of laser drift. If a better frequency fitting were desired additional reference 
etalons could be used. 

Figure H] shows the pump-probe transmission features for the upper hyperfine level 
for each isotope; six sub-Doppler features are clearly seen, with their positions being in 
excellent agreement with the predicted values. 
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Figure 4. Experimentally measured transmission plots for D2 showing saturated- 
absorption/hyperfine pumping spectra of (a) the ^^Rb Fg = 2 line, and (b) the *^Rb 
Fg = 3. The vertical reference lines show the expected peak positions. 

4. Results 

4-1- Effects of hyperfine pumping 

Figure shows a plot for the D2 line of the transmission, T, versus linear detuning, 
A/2tt for a probe intensity of 1.6 /iW/mm^, corresponding to ///sat = 0.1. The expected 
transmission is also plotted. The temperature measured using the thermocouple 
was adjusted at the 0.1°C level in order to fit to the measured data. Reasonable 
agreement is obtained for transitions from the upper hyperfine level of the ground term 
Fg = X+1/2 — > Fe (labelled i and ii), often referred to as the "laser cooling" transitions. 
Poor agreement is seen for transitions from the lower hyperfine level of the ground term 
Fg = I — 1/2 ^ Fe [Hi and iv), often referred to as the "repump" transitions. 

Although the power of the beam is such that the intensity is less than the saturation 
intensity, the assumption that the atomic population has not been influenced by the 
propagation of the probe through the medium is obviously invalid. The process by 
which an atom is excited from one Fg level and is transferred by spontaneous emission 
into the other Fg level is known as optical, or hyperfine, pumping. Allowing for transfer 
out of the two-level system is known to modify the absorption process [ST] . 

Notice for the Doppler-broadened transitions ^^Rb, Fg = 2 ^ Fe (i) and ^^Rb, 
Fg = 3 ^ Fe (ii) in Figure [5] the agreement between theory and experiment is excellent 
on the high-frequency side of the resonance but poor on the low-frequency side - the 
presence of optical pumping not only reduces the peak absorption but also distorts the 
lineshape [32]. 

To investigate this further a sequence of spectra were recorded for different probe 
powers, for both D transitions. Figure [6] shows the line-centre transmission for (a) the 
D2 transition in a room temperature cell, and (b) the Di transition in a cell heated to 
36°C. The laser intensity has been normalized in terms of the saturation intensity |13j . 
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Figure 5. Transmission tlirough tlie vapour cell at 25.4°C, showing the role of 
hyperfine pumping in absorption spectroscopy on the D2 line. The solid red line 
is the experimentally measured transmission and the solid black the predicted value. 
Superimposed are the expected transmissions for the individual hyperfine transitions, 
according to the colour scheme in Figure f . 

Consider the closed hyperfine-resolved transition D2 line: Fg = X+1/2 ^ Fe = X+3/2. 
Owing to the AF = 0, ±1 selection rule atoms excited into this state have to decay to 
the ground state from which they started. These transitions have a significantly larger 
oscillator strength than the two neighbouring transitions Fg = X+1/2 ^ F^ = X±l/2, 
consequently the agreement with the theory which neglects transfer into other ground 
states is good, as was apparent in the high frequency side of the spectrum of the 
transitions i and ii in Figure O For the Doppler-broadened transitions ^^Rb Fg = 2 ^ F^ 
{in) and ^''Rb Fg = 1 ^ Fg (iv) there are two closed transitions, Fg = 1 — 1/2 ^ Fe = 
X — 3/2. However these have similar linestrengths to their neighbouring transitions, and 
hence do not dominate the absorption profile. There are no such closed transitions in 
the Di spectrum. The conclusion therefore is that great care has to be taken to ensure 
that the probe beam intensity is sufficiently low that hyperfine pumping does not occur 
during an atom's transit through the beam - this places a far more strict limitation on 
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Figure 6. (a) The transmission at the centre of the four Doppler-broadened absorption 
features for D2 is pfotted against laser intensity. The data points correspond to the 
measured transmission, and the soHd hues show the transmission expected. The dotted 
hues are guides to the eye. The cell was at 25°C. (b) Similar for Di where six Doppler- 
broadened absorption features are observed. Here the cell was at 36°C in order to 
make the minimum transmission comparable to the D2 line. 



the upper intensity to be used in contrast to the condition / < /gat vahd for two-level 
atoms [S5] . 

An alternative way to visualize the relative importance of optical pumping is to 
plot the normalized absorption coefficient a{I)/a{0). This is done in figure [7] for the 
two D2 repump transitions, and all Di transitions. Also included is the theoretical 
prediction for a Doppler-broadened medium consisting of two-level atoms [13J of the 
form + ///sat- It is apparent that optical pumping reduces the absorption at 

intensities much weaker than those necessary to saturate the transition. The large error- 
bars for small intensities are a consequence of the very low light levels and consequently 
poor signal-to- noise. The data are fit to curves of the form 1 / ^1-1- /5(///sat), with 
(3 being a parameter that characterises the effective reduction in saturation intensity. 
This is done as a guide to the eye, and care should be taken not to over interpret this 
parameterisation. In this work the beams had a fixed width and the intensity was varied 
by changing the probe power. It is possible to realise the same intensity with different 
power beams of different radii; in this case the presence of optical pumping means that 
knowledge of intensity alone is not enough to predict the absorption strength [33] . 

4-2. Comparison of experiment and theory 

With knowledge of how weak the probe beam had to be, we performed a series of 
experiments to test the agreement between our theory for the Doppler-broadened 
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Figure 7. Normalized line-centre absorption showing the effects of hyperfine pumping. 
The solid line is the theoretical prediction for a Doppler-broadened medium of two-level 
atoms. The dotted lines are guides to the eye. Hyperfine pumping on open transitions 
is seen to be more effective at reducing the line-centre absorption than saturation. 

absorption profile of rubidium vapour and experiment. Tlie probe intensity was 
32 nW/mm^, corresponding to ///sat = 0.002. Figure [H] shows transmission spectra 
at tfiree different temperatures (16.5°C, 25.0°C and 36.6°C) for the D2 line. There is 
excellent agreement between theory and experiment; the rms discrepancy is at the 0.2% 
level. Note that the measured absorption is still slightly smaller than the predicted 
value. This could arise due to the broad pedestal of the emission from the laser, and 
also the finite laser linewidth which is of the order of 0.1% of the Doppler width. 

5. Conclusion 

In summary, we have studied Doppler-broadened spectra for the Rb D lines. A model 
was developed which allows the absorption profile and refractive index to be evaluated 
in terms of tabulated functions. Excellent agreement was found between theory and 
experiment for transmissions ranging from 5 to 95%. We showed that the effect of 
hyperfine pumping is significant for open transitions, and outlined how to achieve the 
weak-field limit. The weak-probe limit is only reached for ///sat ~ 0.001 for a beam 
width of 2 mm. Our model allowing quantitative predictions of the absorption and 
dispersion in alkali metal vapour will both aid the burgeoning field of controlled light 
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Figure 8. Transmission plots for the comparison between experiment and theory, at 
temperatures of 16.5°C (top), 25.0°C (middle), and 36.6°C (bottom). Red and black 
lines show measured and expected transmission respectively. Below the main figure is 
a plot of the difference in transmission between theory and experiment for the 16.5°C 
measurement. 



propagation [H [9l [TO] and in the understanding of the spectra obtained in widely used 
laser locking schemes [lEl [TTl [18] . 
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Appendix A. Vapour Pressure and Number Density 

The rubidium vapour cell contains ^^Rb (relative atomic mass Mgs = 84.911789738) 
and ^^Rb (Mgy = 86.909180527) in their natural abundances of 72.17% and 27.83% 
respectively [34j- The vapour pressure (in Torr), p, for solid rubidium is given by the 
following equation [35] . 

logiop = -94.04826 - ^^^^^^^ - 0.03771687 x T + 42.57526 x log^oT, (A.l) 

and for liquid rubidium is given by 
4529 635 

logiop = 15.88253 + 0.00058663 x T - 2.99138 x logioT. (A.2) 
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Using this vapour pressure, the number density, J\f, of rubidium atoms can be calculated, 

(A.3, 

/cbT 

The melting point of rubidium is 39.31 °C. The factor of 133.323 converts the vapour 
pressure from Torr to Pa. Since there are two isotopes present in the cell, the number 
densities need to be calculated separately according to their abundance. 

Appendix B. Transition Coefficients 

The values of the transition strength factors Cp of the D lines are tabulated in this 
Appendix. 



Table Bl. Cj, for the £>2 line of (a) «^Rb, and (b) ^^Rb. 
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Table B2. Cj, for the Di line of (a) ^SRb, and (b) «^Rb. 
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